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We report the results of a search for R parity violating (RPV) interactions leading to the pro- 
duction of supersymmetric sneutrinos decaying into e/x final states using 5.3 fb~ of integrated 
luminosity collected by the DO experiment at the Fermilab Tevatron Collider. Having observed no 
evidence for production of e/x resonances, we set direct bounds on the RPV couplings A311 and A312 
as a function of sneutrino mass. 

PACS numbers: 14.80.Ly, 12.60.Jv,13.85.Rm 



In all standard model (SM) interactions, baryon num- 
ber, B, and lepton number, L, are separately conserved. 
In the supersymmetric (SUSY) extension of the SM, B 
and L violating interactions are generally allowed. A new 
multiplicative quantum number is therefore introduced, 
R = (_i)2S+3(B-z,) ) defined in terms of B, L and the 
spin quantum number S, which distinguishes SM parti- 
cles (R=+l) from their SUSY partners (R = -1) Q. 
If R parity is conserved in the minimal extension of the 
SM, no B and L violating interactions can occur. 
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All SUSY particles are pair-produced if R parity is 
conserved. RPV interactions allow single production of 
SUSY particles which significantly reduces the energy re- 
quired to observe them at a collider. 

The most general renormalizable gauge invariant R 
parity violating (RPV) supersymmetric potential can be 
found in [2J. The terms in the Lagrangian relevant to 
this analysis are 



£rpv 

1 

-A' 



\jk (piLlkRljl — VjLlkRhl) 



i jk \piLdkndji) + h.c, (1) 



where the indices i,j, k — 1,2,3 refer to fermion gen- 
eration; I and d are the SM lepton and down quark 
fields; v is the field of the SUSY partner of the neu- 
trino, the sneutrino. These terms lead to the produc- 



tion of a single sneutrino in dd scattering. The search 
is performed under the hypothesis that only the third 
generation sneutrino {u T ) is produced and that it is the 
lightest SUSY particle. All couplings apart from A 311 
and A312 = A321 = — A231 = — A132 are therefore assumed 
to be zero and the sneutrino decay is determined by the 
e/u and dd modes. 

In this Letter, we report on a search for resonant pro- 
duction of a sneutrino decaying into an electron and a 
muon in pp collisions at y/s = 1.96 TeV at the Fermilab 
Tevatron collider in a data set with an integrated lumi- 
nosity of 5.3 fb -1 , collected between April 2002 and June 
2009. Observation of this process would provide direct 
evidence of physics beyond the SM [3] . Previous searches 
for high-mass e/j, resonances by the CDF [4| and DO [5[ 
Collaborations were based on integrated luminosities of 
1.0 fb -1 . In addition to using a significantly larger data 
sample, the signal sensitivity has been improved by in- 
creasing the lepton acceptance and by applying a neural 
network (NN) selection to distinguish jets from electrons. 

Indirect two standard deviation bounds on the cou- 
pling constants, under the single coupling dominance as- 
sumption with a degenerate sparticle mass spectrum of 
M = Mp T = 100 GeV, are given in Ref. [| as 



A 311 <0.12,A 3 i2<0.07. 



(2) 



The sneutrino production cross section is determined by 
these two couplings and the sneutrino mass M^ T . The 
final state is characterized by an electron and a muon, 
both of which are well-isolated and have high transverse 
momentum, px, of approximately Mi> r /2, and by a peak 
in the invariant e/i mass at M„ T . The dominant SM back- 
ground processes for this event topology are the produc- 
tion of Z/j* -> tt, dibosons (WW/WZ/ZZ), it pairs, 
and W bosons in association with jets, where a jet is 
misidentified as a lepton. 

The DO detector [6| comprises a central tracking 
system in a 2 T superconducting solenoidal magnet, 
surrounded by a central preshower detector (CPS), a 
calorimeter, and a muon system. The tracking system, a 
silicon microstrip tracker (SMT) and a scintillating fiber 
tracker (CFT), provides coverage for charged particles in 
the pseudorapidity range |ry| < 3 Q- The CPS is lo- 
cated immediately before the inner layer of the calorime- 
ter and is formed of approximately one radiation length 
of lead absorber followed by three layers of scintillating 
strips. The calorimeter consists of a central calorimeter 
(CC) covering up to \r/\ f=a 1.1, and two end caps (EC) 
extending coverage to \r/\ « 4.2. Each consists of an in- 
ner electromagnetic (EM) section, followed by a hadronic 
section. The EM calorimeter has four longitudinal layers 
and transverse segmentation of 0.1x0.1 in 77— space, ex- 
cept in the third layer, where it is 0.05 x 0.05. The muon 
system resides beyond the calorimeter and consists of a 
layer of tracking detectors and scintillation trigger coun- 
ters before 1.8 T iron toroidal magnets, followed by two 
similar layers after the toroids. The coverage of the muon 
system is \r]\ < 2. The data acquisition system consists of 



a three-level trigger, designed to accommodate the high 
instantaneous luminosity. For final states containing an 
electron with px > 30 GeV, the trigger efficiency is close 
to 100%. 

To simulate signal kinematics in the DO detector, 
parton level signal events are generated using the 
COMPHEP [8] leading order Monte Carlo generator and 
then processed through pythia [9J to include parton 
showering, hadronization, and particle decays. SM back- 
ground processes are generated with pythia, except for 
the W+jets inclusive samples, which are generated with 
ALPGEN [10| and PYTHIA for parton showering. All sig- 
nal and background simulations use the CTEQ6L1 [111 ] 
parametrization of the parton distribution functions. 

We use next-to-next-to-leading order (NNLO) cross 
sections for Drell-Yan Z/7* -*££(£■ = e , r, u) pro- 
cesses [12| and NLO cross sections for diboson [13| , tt [14J , 
and VF+jets production [15]. All signal and background 
events are processed with a detailed GEANT-based DO de- 
tector simulation [16| and are corrected for trigger effects 
and for the differences in the reconstruction efficiencies 
in the simulation compared to those in data. 

Electrons are selected by requiring an EM cluster in the 
CC or in either EC with trans verse energy Ex > 30 GeV 
within a cone of radius 1Z = ^(A^) 2 + (A77) 2 = 0.2. The 
EM cluster in the CC must be in the range \rj\ < 1.1 and 
in the EC in the range 1.5 < |??| < 3.2. At least 97% 
of the cluster energy must be deposited in the EM sec- 
tion of the calorimeter and the energy must be isolated 
in the calorimeter, [E tot {0A) - £ EM (0.2)]/E: EM (0.2) < 
0.07, where Etot(H) and Eem(R-) are the total energy 
and the energy in the EM section, respectively, within a 
cone of radius 1Z around the electron direction. A track 
must point to the EM cluster for all electron candidates. 
A multi- variable likelihood discriminant, which includes 
information from the spatial track match, must be con- 
sistent with that for an electron in the CC. An NN is 
trained using information from the tracker, the calorime- 
ter, and the CPS to further reject background from jets 
misidentified as electrons. The electron must also be spa- 
tially separated from reconstructed muons. 

For the muon candidate, we require that the associ- 
ated central track pt exceeds 25 GeV and that the time 
measured for hits in the muon scintillation counters is 
consistent with an interaction originating from a pp col- 
lision. The central track fit must have x 2 /ndf < 4 and 
the distance of closest approach (dca) of the track to the 
beam spot in the plane transverse to the beam direction 
should be less than 0.02 cm if the track has SMT hits and 
less than 0.2 cm otherwise. The sum of the transverse 
energy of calorimeter cells in the annulus 0.1 < 1Z < 0.4 
around the muon direction must be less than 2.5 GeV, 
and the sum of the transverse momentum of all tracks 
besides the muon track within 1Z — 0.5 must be less than 
2.5 GeV. 

Events are required to have exactly one high px iso- 
lated electron candidate and one high px isolated muon 
candidate. There is no requirement on the charge of the 
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FIG. 1: (color online). Distributions of (a) $ T and (b) 
A(fi($ T: p) between the direction of the $ T and the muon 
for $ T > 20 GeV. The signal simulation is shown for Mo T = 
100 GeV (400 GeV) and a x BR = 40 fb (12.5 fb). 



two leptons. Furthermore, events are only considered 
if the primary vertex is reconstructed within 60 cm of 
the center of the DO detector in the z coordinate along 
the beam axis and if the difference between the z coordi- 
nates of the muon and the electron at the dca is less than 
2 cm. Background from Z/j* — ¥ tt decays is heavily sup- 
pressed by the high px requirement for the two leptons. 
To reduce the background from tt production, events are 
rejected if they have at least one jet with pt > 25 GeV 
in the range \r)\ < 2.5, where jets are reconstructed using 
an iterative seed-based cone algorithm [17| ■ Figure UJa) 
shows the missing transverse energy ($t) distribution 
and Fig. [T^b) the distribution of azimuthal angles be- 
tween the $ T and the muon direction, A0(j£ T ,/i), for 
events with $ T > 20 GeV. Good agreement between 
data and the total SM predictions is observed. Signal 
events have low p T , but due to the limited muon mo- 
mentum resolution, some $ T is expected in e/i signal 
events that is either pointing in the muon direction or 
opposite to it. This is observed in Fig. QJb) for signal 
events and for the topologically similar Drell-Yan pro- 
cess. Events with $, T > 20 GeV are therefore rejected 
only if 0.7 < A<j>($ T , p) < 2.3. 



FIG. 2: (color online). Distributions of (a) M efi and (b) 
A0(e, p) for data compared to the sum of all background 
processes after the final event selection, rejecting, in addi- 
tion, events with $ T > 20 GeV and 0.7 < A<j)($ T ,p) < 2.3. 
The signal simulation is shown for Mp T = 100 GeV (400 GeV) 
and a x BR = 40 fb (12.5 fb). 



TABLE I: The numbers of selected events in data and esti- 
mated background contributions for an integrated luminosity 
of 5.3 fb -1 with their total uncertainties. 



Process 



Number of events 



Drell-Yan (Zpf) 254 ± 26 

Diboson {WW, WZ, ZZ) 116 ± 12 

it 5.8 ±1.0 

W+jets 34.1 ±5.9 



Total background 



410 ± 38 



Data 



414 



The resulting distribution of the electron and muon 
invariant mass M eft and the azimuthal angle A<^(e, fj.) 
between the electron and muon direction are shown in 
Fig. [5] The kinematic variables of the ep, final state are 
well described by the sum of the SM background contri- 
butions. The background contributions and the number 
of selected candidates are summarized in Table |TJ About 
80% of the Drell-Yan background is due to Z/"f* — » tt 
events. 
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FIG. 3: (color online). The observed and median expected 
upper limits on a x BR for the process pp^f T + X->e^ + X 
as function of Mp T . The median expected limits are shown 
together with the ±1 and ±2 standard deviation bands. The 
theoretical cross sections for A 311 = 0.003 and A312 = 0.005 
and 0.07 are also shown. 



There are 414 candidate events found in the data. The 
expectation from SM processes is 410 ± 38 events, where 
the uncertainty includes the statistical uncertainty of the 
MC, the systematic uncertainties from the integrated lu- 
minosity (6.1%), reconstruction and trigger efficiencies 
(0.5%), which are all taken to be fully correlated be- 
tween the background sources, and the uncertainties on 
the cross sections (Z/j* -> it 3.5%, tt 14.8%, diboson 
production 2.7% - 6.6%, and W+jets 8.5%). Additional 
PDF uncertainties on the signal acceptance are estimated 
from the CTEQ6.1M eigenvector PDF sets and lie in the 
range 0.4% — 0.6%, depending on Mp T . 

The M efl distribution is used to calculate an upper 
limit on the production cross section multiplied by the 
branching ratio, a x BR, for the process pp — > v T + X — > 
e/ji+X with a modified frequentist (CL S ) method 18], un- 
der the assumption that the total width of the produced 
resonance is much narrower than the detector resolution. 
The observed cross section upper limits as a function of 
the My T hypothesis are shown in Fig. [3J together with 
the median expected limits. 

These limits are translated into upper limits on cou- 
plings as a function of Mp T using the theoretical signal 
cross section Q. A mass dependent K- factor, ranging 
from 1.64 at M Vt = 100 GeV to 1.29 at Mp T = 500 GeV, 
is applied to the cross section to include next-to-leading 
order (NLO) QCD corrections [19j . The limits are ob- 
tained by fixing one of the coupling constants and then 



setting the upper limit on the other for different Mp T . 
In Fig. 21 the observed upper limits on A 311 for four as- 
sumed values of A312 are shown. For Mp T — 100 GeV and 
A312 < 0.07, couplings A 311 > 6.2 x 10~ 4 are excluded at 
the 95% C.L. 

In summary, we have searched for a high mass, narrow 
e/it resonance in DO data corresponding to an integrated 
luminosity of 5.3 fb^ 1 . A total of 414 en events is selected 
in data, in agreement with the predicted number of SM 
events, 410 ± 38. The kinematic distributions are well 
described by the SM predictions. The upper limits on the 
production cross section for such a resonance are about a 
factor of 5 lower than for our previous result [5j and about 
a factor 6 lower for M„ T = 500 GeV compared to a search 
performed in the efj, channel by the CDF Collaboration 
with 1 fb _1 of integrated luminosity |4j. We have also 
derived limits on the parameters of a SUSY model that 
predicts a sneutrino resonance decaying into e/i via R 
parity violating production and decay of sneutrinos. 
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FIG. 4: (color online). The 95% C.L. observed upper limits 
on A311 for four values of A312 as a function of Mc T . 
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